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EFFECT  OF  RADIAL  FLOW  BETWEEN  THE  ROTATING  DISC  AND  HOUSING  ON  THEIR 

RESISTANCE  AND  HEAT  TRANSFER 

L.  A,  Dorfraan 


The  problem  under  consideration  was  set  up  In  connection 
with  a  study  of  air  cooling  of  gas  turbine  rotors  by  a  radial 
flow  and  also  of  their  liquid  cooling  by  means  of  shields 
placed  on  the  ends  [-4].  It  Is  necessary  to  Investigate  the 
effect  of  a  flow  on  the  hydrodynamics  of  a  stream  close  to  a 
disc  rotating  In  the  housing  when  studying  the  operation  of 
centrifugal  p\imps,  compressors,  tiirbines,  etc.  [5]. 

Certain  results  of  the  theoretical  Investigation  of  these 
problems  are  cited  below  and  a  comparison  Is  made  with  known 
experimental  data. 


1.  Solution  of  the  basic  equation.  Let  us  examine  the  equation 
of  momentum  for  a  stream  close  to  a  disc  rotating  In  the  housing  (Plg.  l) 
In  the  presence  of  a  radial  flow 

^ |r*  \  v,v^dy^  =  (t,  I,  •  -  T,  Jo)  (  1 ,  i) 

tt 

Here  Vp,  v^  denote  respectively  the  radial  and  peripheral  compon¬ 
ents  of  velocity.  The  magnitudes  of  the  friction  on  the  disc  T^p  \q 
and  on  the  housing  wall  cem  be  presented  In  the  form 

T, !,  =  PP,**,  T,  lo  =  p»o‘*  ( 1  •  2) 
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Let  us  introduce  the  average  discharge  peripheral  velocity  in  the 


gap  v_  ,  so  that 
% 


•  • 

2n  ^  v^v^dy  =  ^  v4y  = 


(1.3) 


Where  Q  is  the  per  second  voltune  discharge  of  the  fluid  in  the  gap. 

Then  Eq.  (l.l)  will  have  the  form 

Q  ^  {r»,o)  =  2nr*  (v*  -  ^  ^ 

The  known  experimental  data  [2,  4]  show  that  for  s'ufficiently 

small  relative  gaps  s:ri  and  for  discharges  through  the  gap  the 

values  of  v„  are  close  to  the  magnitude  of  the  peripheral  velocity 
9o 

v^  in  the  middle  of  the  gap 

1/2  (1,5) 

It  is  necessary  to  note  that 
the  approximate  solution  of  Eq .  (l.l) 
for  a  turbulent  regime  was  first  de¬ 
rived  by  A.  A.  Lomakin  [5].  However, 
in  this  case  the  values  of  the  co¬ 
efficients  of  the  local  resistance 
of  the  disc  and  housing  were  tedcen 
as  constant  along  the  radius  regard- 


e:, 


c* 


Pig.  1.  Diagram  of  the 
problem. 


less  of  the  number  R  =  r^a)/v;  their  numerical  values  were  taken  from 
the  appropriate  experlmantal  data  for  a  disc  in  the  absence  of  a  dis¬ 
charge  flow. 

In  order  to  obtain  a  more  accurate  solution  we  will  use  the  known 
regularities  of  a  tvirbulent  flow  which  were  obtained  while  studying 
streams  in  tubes  and  along  a  flat  wall.  As  is  known  the  use  of  these 
regularities  for  calculating  flows  around  rotating  bodies  yields  re¬ 
sults  which  algree  well  with  experimental  data  [6] . 
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Experiments  with  rotating  discs  In  the  presence  of  a  radial  flow 
In  the  gap  show  [2]  that  there  Is  no  distinctly  outlined  flow  core  and 
we  cannot  speak  about  the  existence  of  separated  boundary  layers  with 
an  Intermediate  potential  core.  Therefore,  we  will  assume  that  start¬ 
ing  with  an  Initial  radius  r^,  the  boundary  layer  at  a  fixed  wall 
closes  with  the  boundary  layer  on  the  disc. 

Taking  Into  consideration  the  comments  made  above  we  can  represent 
the  velocity  profiles  In  the  form  of  an  exponential  law;  close  to  the 
rotating  disc 


(1.6) 


Close  to  the  housing  wall 


p,  _  ^  »,« (» —  y)  > 

p,.  “  1  V 


(1.7) 


Here,  as  Is  known  [7], 

/I  =8.74,  m  = 

Starting  from  the  assumption  (1.5)  we  will  find  the  relation  be- 


teeen  v^^^,  v^^  and  v^^, 


(1.8) 


hence,  we  obtain 


(1.9) 


We  Introduce  the  designation;  z  =  I  m,  x  =  r  I  n  }  after  substitution 
of  (1.9)  Into  (1.4)  we  obtain 

.^(zx*)=  — 2)'^  — *■^1  whM  m  =  l,  >4  =  8.74  (l.lO) 

The  magnitude 

^  =  0.0268  ^  -S’ '  (1*11) 
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Is  denoted  In  terms  of  p. 

Here  Is  the  starting  radius  and  ri  Is  the  final  radius. 

The  value  of  p  shows  on  what  parameters  depends  the  solution  of 
the  problem  ( It  also  depends  on  the  condition  at  the  starting  radius 
X  =  l)  :  when  the  medium  Is  sent  to  the  disc,  as  shown  In  the  diagram 
In  Pig.  1,  z  =  0  when  x  =  Ij  If  the  medium  Is  sent  to  the  gas  along 
the  rotating  shaft,  ^  Is  close  to  unity. 

Prom  Eq.  (l.lO)  we  directly  obtain  the  solutions  for  the  limiting 
cases: 

1)  ■  In  the  absence  of  a  discharge  p  =  «,  so  that 

(1  _  2)7.  =  z\  z  =  0.5 

2)  at  very  large  discharges  0,  then 

zx*  =  const,  z  =  0 

The  result  of  numerical  Integration  of  Eq.  (l.lO)  when 
zo  =  ^  =  0  Is  shown  In  Pig.  2. 

The  approximate  solution  of  Eq.  (l.lO)  can  also  be  obtained  if 
we  note  that 

(1  —  z)yj  —  z’/‘  «  1  —  2z  for  •^0<*.<i0.5 

then  the  solution  of  the  equation  for  an  arbitrary  starting  condition 
z  =  Zo  (when  x  =  l)  will  have  the  form 

z  (z)  =  (X)  +  zo  ^exp  ^  (1  -  x'V.)]  (1.12) 

Where  z*  (x)  Is  the  solution  for  the  starting  condition  z©  =  0, 

z"  (X)  =  -fcexp  \^{i  -  [w^  (*’•*  "  (1.13) 

Calculation  by  this  equation  yields  a  result  close  to  the  calcula¬ 
tion  -by  Eq.  (I.IO)  (Pig.  2). 

Using  Pig.  2  and  Pormula  (1.12)  we  can  find  the  values  of  z  (x) 

_4_ 


for  any  z©. 


Formula  (1.12)  shows  that  the  effect  of  the  Initial  twist  abruptly 
decreases  with  distance  from  the  starting  radius;  this  effect  also 
decreases  with  an  increase  of 


Fig.  2.  The  dependence  of  the  rela¬ 
tive  twist  of  flow  z°  =  V  <p°/rto  in 
the  middle  of  the  gap  on  x  =  r/ro 
for  different  values  of  the  para¬ 
meter  p  at  different  relative  dis¬ 
tances  from  the  axis  of  rotation 
(calculation  by  Formula  (l.lO)); 
the  twist  at  the  starting  radius 
Xi  =  1  is  0  ( zo  =  0,  o  is  calculated 
by  Formula  (1.5)  . 


2.  Resistant  moment.  Com¬ 
parison  of  the  calculations  with 
experimental  data.  Having  deter¬ 
mined  the  values  of  z  ( x)  ,  we 
can  calculate  the  friction  and 
resistant  moment  M  of  the  disc 
surface  being  blown  over 


M 


=  2jt  ^  r*T,  |o  dr 


By  means  of  Formulas  (1.2) 
and  (1.9)  when  m  =  1/7  we  will 
obtain  for  the  coefficient  of 


the  resistant  moment 


CjK  - 


M 


Vipfr 


=  0.337t,-"/.R -/.  _  ^p.dx 


(2.1) 


The  result  of  the  calculations  of  c  when  zo  =  0  is  shown  in 

m 

Fig.  5.  Let  us  compare  the  results  obtained  above  with  the  experimental 
data.  The  calculated  velocity  distributions  on  the  gap  derived  from  ^ 

Formulas  (1.6)  and  (1.7)  are  compared  in  Fig.  4  with  the  experimental  ‘ 

[2].  As  we  see,  as  the  discharge  increases  the  agreement  between  the 
calculations  and  the  experiment  deteriorates.  The  agreement  also 
deteriorates  when  the  relative  gap  Increases.  In  addition  to  the  I 

causes  indicated  in  section  1,  the  nonuniformity  of  the  air  distribution  i 
throughout  the  section  at  the  starting  radius  can  also  have  a  substem-  ^ 
tlal  effect. 
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Pigixre  5  shows  the  calculated  and  experimental  [4]  values  of  the 
relative  magnitude  of  the  peripheral  velocities  In  the  middle  of  the 
gap. 

It  Is  necessary  to  note  that  the  Initial  twist  was  not  measixredj 
the  estimate  based  on  measurements  of  two  relative  radii  x  =  1.3  and 
1.915  yields  a  magnitude  zo  **  O.25. 

When  p  =  00,  l.e..  In  the  absence  of  a  discharge  flow,  the  experi¬ 
mental  values  are  close  to  0.37  Instead  of  0.5  as  calculated, which  Is 
due  to  the  effect  of  the  cylindrical  rim  of  the  housing  and  rim  of 
the  disc  ( see  [6]  )  . 

We  see  from  the  graphs  In  Pig.  5  that  the  effect  of  the  rims  of 
the  housing  and  disc  decreases  with  an  Increase  in  discharge.  Plgure 
5  also  Illustrates  the  fact  that  the  effect  of  the  Initial  twist 
decreases  with  an  increase  In  p  and  with  distance  from  the  starting 
radio . 


Pig.  3.  Dependents  of  the  coefficient  of  the  resistant 
moment  c*  =«„«'/* (Wn)''*  of  the  rotating  disc  on  riiroj 
calculated  by  Pcrmula  ( 21) 

The  Increase  In  the  resistant  moment  of  a  rotating  disc  with  an 
Increase  of  the  discharge  of  the  medium  through  the  gap  was  first 
Investigated  by  V.  S.  Sedach  [2].  Plgure  6  shows  the  Increase  In 
the  resistant  moment  of  a  disc  rotating  In  air;  the  calculated 
values  are  plotted  In  the  same  place.  In  spite  of  the  large  magnitude 
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of  the  relative  gap  there  Is  a  satisfactory  agreement  between  the 
calculation  and  experiment.  The  increase  in  divergence  with  an 
Increase  of  discharge  can  be  explained  by  the  fact  that,  due  to  the 
large  gap  during  the  feed  of  air  through  the  central  tube  the  picture 
of  flow  approaches  the  case  of  an  air  stream  feed  to  the  center  of 
the  rotating  disc. 

0.6 

0.i 

02 

0 

■  ai 


Pig.  4.  Dlstrubutlon  of  peripheral  velocities 
v°  =  V  /  r  03  in  the  gap  between  the  rotating 
disc  and  housing;  the  circles  are  the  experi¬ 
mental  data  of  V.  S.  Sedach  [2],  the  cu3?ves 
are  the  calculated;  in  this  case  Ri  =  10®  , 
s/ri  =  0.055j  a)  K  =130  =  O.OO83) ;  b) 

=  90  =  0.00645^)  . 

The  data  of  V.  S.  Sedach  and  A.  N.  Nespela  [8]  obtained  for  water 
differ  substantlalfer  from  the  calculated  and  experimental  values  obtained 
for  air.  According  to  their  opinion  there  is  an  effect  of  the  mass 
forces  which  are  characterized  by  the  Galilei  number  G  =  grf/v®. 

3.  Heat  transfer  of  the  disc.  To  calculate  the  radial  cooling 
of  a  gas  turbine  disc  rotating  in  a  narrow  chamber  we  must  find  the 
effect  of  the  discharge  on  the  heat-transfer  coefficient  of  the  disc. 

At  a  Prandtl  number  value  of  P  =  1  and  a  quadratic  distribution 
of  the  mean  temperatxire  differences  between  the  disc  and  the  housing, 
as  shown  by  the  author  [6],  there  occurs  a  similarity  of  the  tempera¬ 
ture  distribution  in  the  gap  to  the  distribution  of  the  peripheral 
velocities 
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(7’°  -  T,)  /  (ro  -  f.)  “  * 


(3.1) 


Here  is 
s 

of  the  disc>  T° 


the  temperature  at  a  fixed  wall, 
is  the  temperature  In  the  middle 


Is  the  temperature 
of  the  gap. 


Pig.  5.  Magnitude  of  the  angle  of  the  relative  twist 
of  flow  zo  =  V  /roi  In  the  middle  of  the  gap  on  two 
radii j  1)  calcfllation  at  z©  =  0.25;  2)  calcuatlon  at 
Zo  =  0;  experimental  circles  are  based  on  the  data  of 
B.  P.  Mironov  [4]  ( r©  =  117  mm) 

The  heat  flux  q  from  a  unit  of  disc  surface  equals 


7  —  |o  (fo  —  Tg)  /  no 


(3.2) 


By  means  of  these  formulas  we  can  calculate  the  dimensionless 
coefficient  of  local  heat  transfer  at  radius  r  (Nusselt  number) 


N  =  gr  I  {To-  T°)  X 


(3.3) 


Having  substituted  T^p  |  o  from  (1.2)  Into  (5.2),  we  obtain 


N  =  0.C268  (1  -  (r  /  s)’^* 


(3.4) 


For  the  average  ( Index  m  )  dimensionless  coefficient  of  heat 

transfer  N  we  derive 
m 
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N, 


(To -’n„i 


Here 


r-i 


2« 

n  (rj*  —  ro’) 


T’)  rdr 


Taking  Into  consideration  that  (To-T°)  -v  r®,  we  obtain  finally 


0.0268  (J  -{l  —  zy/‘x'‘/‘dx) 

X 


*1 

^  (1  —  z)x^dx 
i 


(3.5) 


Fig.  6.  Increase  In  the  resistant  moment  AM  kg  cm 
of  a  rotating  disc  at  different  discharges  Q  ks/hr 
of  air  through  the  gap  between  the  disc  and  the 
housing;  1)  calculation  by  Formula  (2.1);  2)  from 
the  experiment  of  V.  S.  Sedach  [2]  (Ri  *  10®,  3/ri= 
*  0.1,  ri  =  0.2  m,  ri/ro  =  4) 


Pig.  7.  Heat  transfer  of  disc  N*  =  (s/rO^/^' 

rotating  In  the  housing  at  P  =  1  and  a  quadratic  distri¬ 
bution  of  the  mean  temperature  differences  (calculation 
by  Formula  (3*5)) 
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Figure  7  shows  the  result  of  rough  calcuatlon  by  this  formula. 

Heat  transfer  of  discs  of  gas  turbines  with  radial  cooling  was 
Investigated  by  V.  M.  Kapinos  [3].  Although  he  measured  the  heat 
transfer  of  the  disc  not  by  calorimetry  and  the  heat  transfer  coeffi¬ 
cient  found  are  in  a  certain  sense  conditional,  it  is  nevertheless 
advantages  to  compare  them  with  the  calculation.  In  the  experiments 
the  distribution  of  the  mean  temperature  differences  had  a  qiiadratlc 
character,  therefore  we  can  use  directly  Fig.  7  if  we  still  take  into 
account  the  effect  of  P  number.  This  effect  can  be  determined  approxi¬ 
mately  from  the  data  for  a  free  disc  [6]  with  the  help  of  the  multiplier 
P°‘®. 


For  the  A  series  of  experiments  a  good  agreement  with  the  calcu¬ 
lated  value  is  obtained  (Pig.  8).  On  doubling  the  axial  gap  (A^  series 
with  a  thin  disc)  the  experimental  values  substantially  exceed  those 
calculated  which,  in  addition  to  making  the  air  distribution  wors^ 
is  associated  mainly  with  the  conditionality  of  the  experimental 


The  experiments  confirmed  the 
calculated  character  of  dependence 

on  Xi  =  rj/ro  and  on  K^.  For 
example  a  change  of  xi  from  2  to 
4  (for  P  =  0.2)  yields,  according 


Pig.  8.  Comparison  of  the 
heat  coefficient  of  a  rotating 
disc;  1])  based  on  the  experi¬ 
ments  of  V.M.  Kapinos  (series 
a)  ;  2)  based  on  formula  (3*5) 
(ti  »  0.32  m,  s  =  10mm, 

=  2,  xi  =  3.8) 


to  Pig.  7,  a  drop  by  24.5^,  and 

according  to  the  experiments 
Xi“°*  ®)  decreases  by  23^. 
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Further  when  Ri  =  10®,  ri/s  =  J>2,  Xi  »  Z>.Qi  a  fourfold  Increase  in 
(from  2to  8)  yields^  according  to  the  calculation^  a  decrease  of 
(by  8^)  and  a  decrease  (by  \K%)  Is  also  obtained  from  the  experiments. 

Therefore,  the  calculated  values  are  In  complete  agreement  with 
the  experimental.  For  a  more  accurate  determination  or  the  agreement 
between  them,  we  must  measure  carefully  the  heat  transfer  coefficient 
of  the  disc  by  direct  calorimetry. 


4.  Heat  reception  and  Shields. 
The  heat  flux  q_  at  the  housing 

D 

wall  when  P  =  1  and  a  qxiadratlc 
distribution  of  the  mean  tempera¬ 
ture  differences  along  the  radius 
Is  determined  in  the  same  manner 
as  that  done  to  the  disc: 


Fig.  9.  A  comparison  of  the  „  ri: 

calculated  N  and  the  experimental  ^$)lra>  I*:. 

N°  of  the  local  heat  transfer  coef¬ 
ficient;  calculation  is  based  on 
formula  (4.20)  and  the  experiments 
are  based  on  the  data  of  B.  P. 

Mironov  for  two  radii  r  =  0.206  m 
and  r  =  0.245  m  ( ro  =117  m  ) 

Then  taking  Into  account  (5*l)  and  by  means  of  (1.2)  and  (1.9) 
we  obtain 


N.  -  0.0268  (r / s)’'*  ^ 

The  extensive  data  on  the  heat  reception  of  shields  in  the  presence 
of  a  discharge  flizx  were  obtained  by  B.  P.  Mlrono  [4] .  A  comparison 
of  the  experimental  and  a  calculated  values  of  local  coefficients  N 
is  made  in  Pig.  9.  For  this  p\irpose  the  calculation  was  performed 
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by  Formula  (4.2)  without  consideration  of  the  effect  of  the  P  number 
and  the  radiation  exchange  was  not  excluded  from  the  experimental 
value.  Because  the  effect  of  the  P  number  can  be  taken  into  account 
approximately  by  the  multiplier  P°'®  =  0.7°’®  =  O.806  and  the  radiant 
heat  based  on  the  experiments  was  about  10^,  we  can  consider  that  there 
Is  a  satisfactory  agreement  between  the  experimental  and  calculated 
data.  On  an  Increase  of  the  gap,  as  Is  apparent,  the  calculated  values 
begin  to  lag  behind  the  experimental.  The  same  result  is  obtained 
with  an  Increase  of  the  discharge,  partlcularily  when  <20. 

In  conclusion  we  will  cite  another  formula  for  the  average  coef¬ 
ficient  of  the  heat  reception  of  the  shield;  It  has  a  form  similar 
to  Formula  (3.5) 

1 

We  note  that  analogous  calculations  can  be  performed  In  the  presence 
of  a  discharge  flow  In  the  gap  between  the  rotating  disc  and  housing 
from  the  peripheral  radius  toward  the  center  of  the  disc. 
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